Formation of molecules near a Feshbach resonance in a ID optical lattice 
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We calculate the binding energy of two atoms interacting near a Feshbach resonance in the 
presence of a ID periodic potential. The critical value of the scattering length needed to produce 
a molecule as well as the value of the molecular binding energy in the unitarity limit of infinite 
scattering length are calculated as a function of the intensity of the laser field generating the periodic 
potential. The Bloch bandwidth and the effective mass of molecules are shown to depend strongly 
on the value of the scattering length due to the correlated motion of the two atoms. 
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The availability of tunable periodic potentials of opti- 
cal nature (see, for example, Ref. 1]) and the quantum 
manipulation of the interatomic forces via Feshbach res- 
onances are opening new perspectives in the field of 
ultracold atomic gases. In the absence of the external 
potential molecules can be formed near resonance only 
for positive values of the s-wave scattering length a, the 
corresponding binding energy being given by the formula 
Eb = —h 2 /ma 2 where m is the mass of each atom. These 
molecules have already been observed experimentally by 
several groups 0, lil HI El ■ An interesting and very rich 
problem is the formation of molecules in an optical lat- 
tice generated by laser fields. Recently, Fedichev et al. 
investigated the case of a 3D tight lattice and showed 
that the bound state can be formed also for negative scat- 
tering lengths down to a critical value. Their approach 
is valid provided that the binding energy is small com- 
pared to the width of the lowest Bloch band so that the 
molecular wavefunction extends over many lattice sites. 

In this Letter we discuss the formation of molecules in 
a ID optical lattice for arbitrary values of the scattering 
length and the laser intensity. If the barriers separating 
two consecutive wells are very high (high laser intensity) 
the system behaves like a scries of 2D discs which lie 
at the bottom of the wells and feel a harmonic poten- 
tial in the direction of the laser field. In this asymptotic 
regime the two-body problem is significantly simplified 
by the possibility of separating the center of mass and 
the relative coordinates of the two interacting atoms ||. 
When this condition is violated many interesting ques- 
tions can be addressed: what is the critical value of the 
scattering length required to ensure the formation of a 
molecule? What is the value of the binding energy of the 
molecule in the so called unitarity limit where the scat- 
tering length approaches an infinite value? How do the 
tunneling properties (bandwidth and effective mass) of 
molecules depend on the value of the scattering length? 
The purpose of this Letter is to provide an explicit answer 
to these questions. 



Let us start by writing the 3D Schrodinger equation 
for two atoms interacting in the presence of a ID periodic 
potential V op t(z): 
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where V(Z,z) — V op t{zi) + V op t{z2) and we have in- 
troduced the center of mass and relative coordinates 
Z = (zi + z 2 )/2 and r = ri - r 2 . 

In Eq. Q the interaction is modeled via a s-wave pseu- 
dopotential with coupling constant g — 4irh 2 a/m, where 
a is the scattering length in vacuum. The solution of 
Eq.JD can be written, for g ^ 0, as 

*(r,Z) = J dZ'G E (r,Z;O,Z')g-£j{T>*{r',Z')) r ,^ , 

(2) 

where G E is the free (g — 0) Green's function associ- 
ated to Eq.QJ. At short distance, the Green's function 
is determined only by the kinetic energy. Neglecting the 
external potential and setting E = in Eq.QJ, we find 
G E=0 (r, Z; 0, Z') = -{m/^ 2 h 2 )[{Z - Z') 2 + (r/2) 2 ]- 1 . 
Therefore the Green's function has the following expan- 
sion 



G E (r, Z;0,Z') 
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-5(Z- Z')+K E (Z, Z') + 0(r), 

(3) 

where the regular kernel K E (Z, Z') depends on the exter- 
nal potential and we used the representation of the delta 
function S(x) = lim r ^o r/ir(r 2 +x 2 ). When inserted into 
Eq.(0l, Eq.@ implies, for r — > 0, the boundary condition 
*(r, Z) — (1/r- l/a)f(Z), where f(Z) is a function of 
the center of mass position satisfying the integral equa- 
tion 



-f(Z) = f dZ' K E (Z,Z')f(Z'). 
9 J 



(4) 



If the center of mass and the relative motion decouple 
the function f(Z) must be an eigenstate of the center of 



mass Hamiltonian independent of the scattering length. 
By inserting it into Eq.©, the dependence on Z factors 
out and we are left with an algebraic equation yielding 
the energy as a function of the scattering length. For 
instance, in the absence of the optical lattice the kernel 
is given by 

K%{Z,Z') = £^J ^-e^ z - z 'W\E\ + VPy4m. 

(5) 

The lowest energy solution of Eq.© corresponds to f(Z) 
constant and yields E = — ft 2 /ma 2 . 

For periodic potentials, the kernel has the property 
Ke(Z, Z') = Ke(Z + jd, Z' + jd) where j is any integer 
and d is the lattice spacing. This means that the general 
solution of Eq.© has the Bloch form f(Z) = e lQZ f Q (Z), 
where /q{Z) is a periodic function, Q being the quasi- 
momentum of the molecule. In the following we shall 
restrict to the lowest bound state given by Eq.©. The 
energy as a function of the scattering length is then fixed 
by the condition that the lowest eigenvalue of the kernel 
K E (Z, Z') in Eq.© is equal to l/g. 

Due to the singular term in Eq. © , the numerical eval- 
uation of the kernel Ke(Z, Z') from the Green's function 
Ge is not trivial. The latter can be written in the general 
form 

4>n 1 ,q 1 {Z)4>n 2 ,q2 (Z)4'n 1 ,q 1 (Z')<j>n 2 ,q 2 

E - e ni {qi) - e n2 {q 2 ) -ft 2 k\/m 

where qs = ft/d is the Bragg wave-vector and <f> n ,q{zi) 
are the eigenstates of H z = — {ft 2 /2m)d 2 /dz 2 + V op t(zi) 
with energies e n (q)- In the following the periodic poten- 
tial will be taken of the form 

V opt { Zl ) = sE R sm 2 (^i), (7) 

as usually provided by two counterpropagating laser 
beams (optical lattices), where s is the laser inten- 
sity, the spacing d being related to the wavelength of 
the laser fields and Er = ft 2 ir 2 /2md 2 is the recoil en- 
ergy. Since the singular term in the rhs of Eq.© 
does not depend on the external potential, we add and 
substract from Ge, the Green function © evaluated 
for s = 0. This permits us to write the kernel as 
K E (Z, Z>) = lim^ [G £ (r, Z; 0, Z') - G° E (r, Z; 0, Z')\ + 
K%(Z,Z'), with K E defined in Eq.©. After these 
manipulations, the integration over k_L in Eq.© con- 
verges and can be performed analytically. The remain- 
ing summation over band indices and the integration 
over quasi-momentum are done numerically. We define 
the binding energy as Eb{Q) — E{Q) — E re f(Q), where 
Eref(Q) — 2ei(Q/2) is the lowest energy state for two 
non interacting particles with total quasi-momentum Q 
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FIG. 1: Binding energy versus inverse scattering length for 
different values of the laser intensity: from top to bottom 
s = 20,10,5,0 (solid line). Also shown are the asymptotic 
behaviour for large d/a (dotted line) and the binding energy 
in harmonic approximation (dashed line) for s — 20. 
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FIG. 2: Critical value of the inverse scattering length as a 
function of the laser intensity. 



and ei is the dispersion of the lowest Bloch band. With 
this notation, E^ is always negative for a bound state. 

The binding energy for quasimomentum Q = as a 
function of the inverse scattering length is shown in Fig^ 
for different values of the laser intensity. 

When the binding energy is large with respect to 
the depth of the lattice (ft 2 /ma 2 ^> sEr), the peri- 
odic potential Q can be treated as a perturbation and 
the wavefunction takes the form *£(r,Z) — (p(Z)<f>b(r) , 
with 4>b(r) = (\ / \f2~Tia)e~ r I a / r . Within first order per- 
turbation theory, the effective potential acting on the 
center of mass of the molecule is given by U e ff(Z) = 
f$(r)V{Z lZ )d 3 r, yielding 

U eff {Z) = sE R - sE R —aTctan(^)cos(2q B Z). (8) 
aqs 2 

We see from Eq.© that the depth of the effective poten- 
tial depends on the scattering length. In the limit a <C d, 
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FIG. 3: Binding energy at unitarity (d/a = 0) as a function of 
the laser intensity (solid line). Also shown are the harmonic 
approximation (dashed line) and the inclusion of anharmonic 
corrections (solid circles). 



the external potential is constant on the scale of the size 
of the molecule and one finds U e ff(Z) = 2V op t{Z). When 
a increases, the size of the molecule becomes comparable 
to the lattice spacing and the effective potential © is re- 
duced. For h 2 /ma 2 3> sE Bl the binding energy is given 
by 



E b (Q) = -h 2 /ma 2 + E'(Q) - E ref (Q), 



(9) 



where E' can be found by solving the Schrodinger equa- 
tion 



— V 2 z + U eff (Z) v(Z) = E' i p(Z) 



4m 



(10) 



The asymptotic behaviour based in Eqs JSJ-JTHJl is plot- 
ted in Fig^with dotted line for different value of s. 

When the laser intensity becomes large (s> 1) and the 
binding energy is large compared to the width of the low- 
est Bloch band, the system enters the quasi-2D regime. 
In this limit, the two interacting atoms are localized at 
the bottom of the same optical well where, to a first ap- 
proximation, the potential Q is harmonic with frequency 
luq = 2y/sEji/H. In this case, the center of mass and the 
relative motion decouple and the scattering problem can 
be solved analytically @, • The corresponding result for 
the binding energy as a function of the inverse scattering 
length is shown in Fig^for s — 20 (dashed line). We see 
that the harmonic approximation slightly overestimates 
the correct binding energy. 

An important consequence of the presence of the opti- 
cal lattice is the formation of a bound state also for nega- 
tive scattering length. However, unlike in the asymptotic 
2D limit, the bound state disappears for values of the in- 
verse scattering length below a critical value d/a cr . The 
dependence of d/a cr as a function of the laser intensity 
is shown in Fig|3 We see that d/a cr = for s = and 
d/a cr — > — oo as s — > oo. 



The behaviour of the binding energy as a function of 
the scattering length sufficiently close to the critical point 
is of the form 1/a — l/a cr = \E b \m* / hC , where m* 
is the atomic effective mass evaluated at the bottom of 
the lowest Bloch band and C — C(s) is a dimensionless 
function of the laser intensity. For s — one has C = 1 
while in tight binding limit one finds C = d/^/2na, where 
a is the width of the variational gaussian ansatz w(z) = 
cxp(— z 2 /2cr 2 )/7r 1 / 4 cr 1 / 2 for the Bloch function 4>i, q (z) ~ 
J2e e tqld w(z~ Id). This is just the renormalization factor 
of the coupling constant g obtained by integration over 
the gaussian profile [ljj. The above results show that 
the optical lattice gives rise to an effective shift of the 
resonance, in analogy with the case of a 3D tight optical 
lattice investigated in Ref.0- 

It is interesting to discuss explicitly the behaviour of 
the binding energy at resonance (1/a = 0) as a function 
of the laser intensity. The numerical result is shown in 
FigED (solid line) . For large s the harmonic approxima- 
tion gives E£°(l/a = 0) = -0.2Uhuj Q = -0A88^/sE R 
(dashed line). The remaining shift AE in the binding is 
due to the non harmonic terms. We have found that the 
inclusion of such corrections via first order perturbation 
theory results in an s-independent shift AE = 0.198i?R 
which well reproduces the full numerical result for s > 10 
(see dots in Fig|3J). 

So far we have discussed the results for the binding 
energy corresponding to vanishing values of the quasi- 
momentum. Differently from the case of harmonic trap- 
ping, where the binding energy does not depend on the 
center of mass motion, in the presence of an optical lat- 
tice the binding energy of the molecule depends on its 
quasi-momentum. This is plotted in Fig^l for s = 2.5 
and different value of the scattering length. We see that 
the modulus of the binding energy increases when Q in- 
creases. Moreover, when the scattering length is nega- 
tive and d/a crosses the critical value d/a cr , the bound 
state with zero quasi-momentum breaks first. The 
bandwidth, defined as w — E(qB) — E(0), is plotted in 
FigGD (Left Panel) as a function of d/a for s = 2.5. For 
a given value of the laser intensity, the curve stops at 
d/a = d/a cr where w = 2e 1 (q B /2) + E b (q B ) < 2t 1 (q B /2) 
because the binding energy at Q = q B is finite at the 
critical point (see FigBJ. The bandwidth for molecules 
is always smaller than the one for single atoms. For in- 
stance for s = 2.5 the width of the lowest band for atoms 
is 0.52Er. This effect results from the correlated mo- 
tion of the two constituent atoms. This is best seen for 
h 2 /ma 2 ^> sEr, where the depth of the effective poten- 
tial felt by the molecule increases as d/a increases. In 
the limit d/a > 1, U e ff(Z) = As E Rm cos 2 (nZ/d), where 
ER m = Er/2 is the recoil energy for a particle with 
mass 2m, yielding an effective laser intensity s m = 4s. 
In particular, in tight binding approximation the atomic 
bandwidth is proportional to the WKB tunneling expo- 
nent e _s °, where So = 7r 2 s 1 ^ 2 /4, while the bandwidth for 
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FIG. 4: Binding energy dispersion as a function of quasi- 
momentum Q for fixed s = 2.5 and for different values of 
the inverse scattering length: from top to bottom d/a = 
2,0, — 1.17(= d/a C r), — 1-66. The solid squares give the dis- 
persion at d/a — 2 evaluated from Eqs IHt- HlUI . 
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FIG. 5: Bandwidth (Left Panel) and effective mass ratio 
2m/ M* (Right Panel) of the molecule as a function of the 
inverse scattering length for s = 2.5 (solid line). The band- 
width at the critical point is always smaller than 2e\(qs/2) 
(dot-dashed line) [See text]. Conversely, at the critical point 
the ratio 2m/ M* coincides with the effective mass ratio m/m* 
for single atoms (dot-dashed line) . The asymptotic behaviour 
for large positive d/a is also shown (dotted line). 



molecules is much smaller, being proportional to e 2S ° . 

In Fig [SI (Right Panel) we show the results for the ef- 
fective mass of the molecule, defined through 1/M* = 
(d 2 E(Q) /dQ 2 )\Q =0 , as a function of d/a for laser inten- 
sity s = 2.5. We see that, in complete analogy with the 
bandwidth, the ratio 2m/M* decreases as d/a increases. 
At the critical point where the molecule breaks, one finds 
M* = 2m*, where m* is the effective mass of a single 
atom. 

In the last part of the Letter we discuss the behaviour 
of the size of the molecule as a function of the scattering 
length and the laser intensity. The typical size of the 



molecule in the radial direction is fixed by the binding 
energy according to Aj^ w h/ y/\Eb\m, where Eb is fixed 
by the scattering length [see Fig.JU]. Due to the optical 
confinement, the molecule is more squeezed in the ax- 
ial direction, so generally X z < X±. For large d/a the 
molecule is not distorted by the optical lattice and there- 
fore Ai ~ A z « a. When d/a decreases, the wavefunction 
becomes more elongated in the radial direction. As long 
as the binding energy is large compared to the width of 
the lowest Bloch band, the two atoms stay in the same 
optical well and therefore X z < d. When this condition 
is violated, the two atoms can hop to different lattice 
sites and the wavefunction becomes spread also in the 
z direction. Finally, close to the critical point, one has 
X±_ i X z 3> d with the anisotropy ratio fixed by the simple 
relation X±/X z — y/m* /m. 

In conclusion, we have presented an exact numerical 
method to calculate the binding energy and the tunnel- 
ing properties of a molecule formed by two atoms inter- 
acting via a short range potential, in the presence of a 
1 D optical lattice, for arbitrary values of laser intensity, 
scattering length and quasi-momentum of the molecule. 
The theoretical predictions could be checked in experi- 
ments on dilute ultracold gases. They could also have 
important consequences on the many-body properties of 
these systems. In particular the formation of molecules 
induced by the periodic optical field is expected to modify 
the behaviour of the BEC-BCS crossover in a degenerate 
Fermi gas. 
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